Tunneling in suspended carbon nanotubes assisted by longitudinal phonons 
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Current-voltage characteristics of suspended single-wall carbon nanotube quantum dots show a 
series of steps equally spaced in voltage. The energy scale of this harmonic, low-energy excitation 
spectrum is consistent with that of the longitudinal \ow-k phonon mode (stretching mode) in the 
nanotube. Agreement is found with a Franck- Condon-based model in which the phonon-assisted 
tunneling process is modeled as a coupling of electronic levels to underdamped quantum harmonic 
oscillators. Comparison with this model indicates a rather strong electron-phonon coupling factor 
of order unity. 

PACS numbers: 



In nano-electromechanical systems (NEMS), mechan- 
ical motion affects electrical current and vice versa [J, 
HI d Of special interest is the study of electron- 
phonon coupling in these devices since tunneling of a 
single electron mav induce a displacement of the mov- 
able structure [3, |y, til lj 111 LLi| • The interaction strength 
is characterized by the dimensionless electron-phonon (e- 
ph) coupling constant g, which is proportional to the ra- 
tio of the classical and the quantum displacement. In 
bulk systems the e-ph coupling is generally weak and the 
coupling constant is orders of magnitude smaller than 
one. However, since the coupling dramatically increases 
with decreasing device mass, NEM-devices may exhibit 
an intermediate to strong e-ph coupling [lj], LLJ, []jj . In 
this regime, current-voltage characteristics are expected 
to exhibit additional steps whose height can be used as 
an estimate of g. For example, g is around one in the C6o 
molecular devices of Ref. [lj| , while measurements on dif- 
ferent C140 samples [l^ indicate a value of g between 0.2 
and 8. 

Carbon nanotubes (NTs) are ideal systems for explor- 
ing electro-mechanical effects since they have a small di- 
ameter, a low mass, and can be defect free on a molecular 
level. In experiments on suspended nanotubes, differ- 
ent methods have already been used to probe the bend- 
ing [30 and radial breathing mode (RBM) 16]. The 
measurements show that the free-hanging tubes operate 
in the underdamped regime of low dissipation. For the 
fundamental bending mode the reported quality factor is 
about 100; for the RBM it is estimated as high as 10000. 

In this Letter we present electronic transport spec- 
troscopy measurements on suspended single-wall nan- 
otubes, which show signatures of phonon-assisted tunnel- 
ing, evidenced by the presence of a series of steps in the 
I — V characteristics. Such steps form a harmonic low- 
energy spectrum, whose energy scale and length depen- 
dence are consistent with that of the longitudinal stretch- 
ing mode. Comparison with the Franck- Condon theory 
shows that the e-ph coupling constant is of order one. 

Devices are fabricated by locating individual nan- 




FIG. 1: a) Schematic drawing of a suspended nanotube 
clamped between two Cr/Au electrodes on top of silicon ox- 
ide. The underlying oxide is partially removed by a wet etch 
step leaving the nanotube suspended. The highly doped sil- 
icon plane is used as a global gate to tune the electrostatic 
potential of the nanotube. b) Scanning-electron microscope 
micrograph of a suspended nanotube. The scale-bar repre- 
sents 200 nm. 



otubes (laser ablation and CVD) on a Si/Si02 substrate 
using an atomic force microscope (AFM) with respect 
to predefined markers. Subsequently, the electrodes are 
made using conventional e-beam lithography techniques 
and thermal evaporation of Cr (5 nm) and Au (50 nm). 
The nanotubes are suspended by removing the under- 
lying Si02 in a wet etch step using buffered HF [~L7^ . 
A schematic sample geometry and SEM micrograph are 
shown in Fig. 1. In the experiment the source and gate 
voltage are denned with respect to the drain, which is 
connected to ground. 

In Fig. 2 we show stability diagrams for three nan- 
otubes measured at 10 mK (a) and 300 mK (b,c) where 
the differential conductance, dl/dV, is plotted versus 
bias and gate voltage. The three metallic nanotubes have 
a length between source and drain contacts, L, ranging 
from 0.14 to 1.2 /im. Their diameter d is between 1 and 
1.4 nm as determined from AFM imaging. In the dia- 
mond shaped regions (Coulomb diamonds) the current is 
zero due to Coulomb blockade, and the charge number 
in the dot is fixed. Regular and closing Coulomb dia- 
monds indicate single dot behavior [lMllf^ in all three 
samples for the gate range shown. Notice that the dia- 
monds in Fig. 2a do close, as shown in the inset, which 
was taken at a higher temperature (300 mK) in a different 
cooldown. The low-bias current, however, is suppressed 
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FIG. 2: Stability diagrams for three different suspended nan- 
otubes with a length in between contacts of 1.2 /mi, 420 nm, 
and 140 nm for (a), (b), and (c) respectively. The conduc- 
tance (dl/dV) is plotted as a function of source-drain volt- 
age, V, and gate voltage, Vg- Blue corresponds to low and 
red to high conductance. Measurements have been performed 
at T=300 mK except in (a), where the base temperature was 
10 mK. (a) Small region of a stability diagram showing closely 
spaced sets of lines running parallel to the Coulomb diamond 
edges for two charge states. At low bias, a strong suppres- 
sion of the conductance is present. Red lines indicate positive 
differential conductance; blue lines negative differential con- 
ductance. Inset: regular diamonds that close are observed in 
a different cool down at T— 300 mK. (b) and (c) Diamond 
crossings for two other samples, again showing lines parallel 
to the diamond edges with energy separations smaller than 
expected for electronic excitations. 



which could be a signature of strong electron-phonon cou- 
pling i 10- 

Excitations of a quantum dot appear as lines running 
parallel to the Coulomb diamond edges in the stability 
diagrams 18]. At such a line, a new electronic level be- 
comes resonant with the leads and an additional trans- 
port channel opens up. The energy of an excitation can 
be determined by reading off the intersection point be- 
tween the excitation line and the Coulomb diamond edge 
on the bias axis |l9|- Furthermore, the excitations corre- 
spond to the charge state of the Coulomb diamond they 
end up in. Electronic excitations in nanotubes typically 
differ between adjacent charge states |20j. In Fig. 2a, 
a dense set of equally spaced excitation lines (starting 
from the first electronic excitation) is clearly visible near 
Vg = 210 mV and 230 mV, i.e., adjacent charge states 
exhibit a similar set of excitations with approximately 
the same energy spacing. The fact that excitations occur 
primarily in one direction is due to asymmetric tunnel 
barriers |l9|- 

The energy differences between the excitation lines of 
Fig. 2 are shown in the insets of Fig. 3. In all three 



cases, the excitation energy is an integer multiple of the 
first (fundamental) excitation. Thus, they form a har- 
monic spectrum with up to 5 levels. A linear fit yields 
an excitation energy of 140, 690, and 530 /ieV for the 
tubes with length 1.2 pm (a), 420 nm (b), and 140 nm 
(c), respectively. These values are an order of magnitude 
smaller than the expected mean electronic level spacing 
given by A = hvp/^L with h the Planck constant and 
vf = 8.1 • 10 5 m/s |2l| the Fermi velocity. 

A more natural explanation for the observed harmonic 
spectra is a vibrational mode coupled to electron tun- 
neling [ll|. Multiple steps with identical spacing would 
than arise from the excitation of an integer number of 
vibrational quanta. Indeed, the observed equidistant en- 
ergy separation is consistent with that expected from the 
longitudinal stretching mode in the nanotubes. In Fig. 4, 
we plot the energy of important low-energy vibrational 
modes of single- wall nanotubes |22|, |23| . For comparison, 
we plot the mean electronic energy level separation, A, in 
black. Blue squares correspond to the fundamental vibra- 
tional excitation energy extracted from the linear fits in 
the insets of Fig. 3. The energy of the radial breathing 
mode (green) does not depend on the nanotube length 
and equals 28 me V/ d( nm) . The bending mode (red) has 
a L~ 2 dependence |25[, and an energy much smaller than 
the measured excitation energy. The stretching mode 
vibration energy (blue) is inversely proportional to the 
length |24j. E = (nh/L)y/Y/p m , where Y is Young's 
modulus, p m is the density and n is the vibrational quan- 
tum number. For nanotubes with p m = 1.3 g/cm 3 , 
Y = 1 TPa the vibrational energy corresponding to the 
fundamental mode is ~ 110 peV/L(pm) |2j. As Fig. 4 
shows, the data are in good agreement with these pre- 
dicted values. 

The coupling of electronic levels with vibrational 
modes (quantum harmonic oscillators) can be described 
in terms of the Franck-Condon model [2^. According to 
the Franck-Condon principle, an electron in an electronic 
transition moves so fast that the nuclear positions are vir- 
tually the same immediately before and after the tran- 
sition. As a consequence, the transition rate is propor- 
tional to the Franck-Condon factors denned as the square 
of the overlap integral between the vibrational wavefunc- 
tions of the two states involved. An important parameter 
is the electron-phonon coupling factor, g = \{^) 2 ''• This 
is the ratio of the classical displacement length, x, to the 
quantum mechanical oscillator length, xq = y/h/mco. Al- 



ternatively, g 



where F is the force on, m the 
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mass of, and uo the frequency of the oscillator. 

For low damping, the vibrational levels remain sharp 
and the Franck-Condon model predicts steps in the 
current-voltage characteristics, that are equally spaced 
in energy (bias voltage). In the presence of strong re- 
laxation, the normalized step heights are given by pj: 
Pn — e~ 9 g n /n\. In the strong coupling (g ^> 1) limit, 
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FIG. 3: Current as a function of source-drain voltage at a voltage indicated by the green lines in Fig. 2. The red lines represent 
the step heights calculated in the Franck-Condon model (see text) for an electron-phonon coupling parameter of 0.95, 1.1, and 
0.5 for (a), (b), and (c) respectively. In the insets, the energy separation between the peaks or steps (lines in Fig. 2) is plotted, 
showing equally spaced, harmonic spectrum. The slope of the drawn line is 140, 690, and 530 fieV for the insets of (a), (b), 
and (c) respectively. 
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FIG. 4: Energy scales of different vibrations and electronic 
excitations plotted on a log scale for a nanotube with a 1.4 nm 
diameter. The radial breathing mode (green) does not depend 
on the length L. The bending mode vibrations (red) have a 
L~ 2 dependence. The mean electronic level spacing (black) 
and the stretching mode (blue) vibrations depend inversely 
on the length. 



the height of the first steps is exponentially suppressed 
(phonon blockade) 0, 0, . Multiple steps only arise 
if g is of the order of one or larger and the observation 
of a spectrum of equally spaced excitation lines there- 
fore indicates that the e-ph coupling in our suspended 
nanotubes must be rather strong. 

In Fig. 3, the red curves represent the step heights (P n ) 
given by the Franck-Condon model with strong relax- 
ation discussed above. The symbols are the experimental 
curves taken at the green lines in Fig. 2. Considering the 
simplicity of the model, reasonable agreement is obtained 
in all three cases. The comparison yields an estimate of 
g of 0.95, 1.1, and 0.5 in Fig. 3a, b, and c respectively, 
indicating that it is approximately length independent. 
We have also performed a similar analysis at other gate 
voltages yielding the same g- values. 

The theoretical curves in Fig. 3 do not exactly fol- 
low the measured ones. Better fits may be obtained 



if the influence of a gate voltage and asymmetric cou- 
pling is considered p or if coupling to excited electronic 
states [22I is considered or if the influence of damping or 
non-equilibrium phonons (weak relaxation) is taken into 
account. In the latter case the peak heights are expected 
to display a non-systematic dependence on g and peak 
number p. Consideration of these effects is, however, 
beyond the scope of this paper. 

The high value of the e-ph coupling (g ~ 1) is remark- 
ably since in graphite the coupling between electrons and 
longitudinal phonons is weak. A source of intermediate 
to strong coupling could be the interaction between longi- 
tudinal and transverse vibrations |23|. Alternatively, we 
find that in suspended nanotubes the same e-ph coupling 
mechanism as in the bulk 29] can lead to a g ~ 1 if the 
electron density is inhomogeneous. The calculation pro- 
ceeds as follows: The interaction energy of electrons with 
the polarization charge is characterized by the energy 
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Here, p(x) is the density of excess charge produced by 
one electron, K(x-x') is an interaction kernel, which we 
approximate by S(x — x') for the case that interactions 
are effectively screened by the gate, and P(x) « epoz(x) 
is the polarization vector. The quantity po ~ hp is 
the total electron density and z(x) is the displacement, 
which in the single- mode approximation becomes z(x) — 
A n sm(7rnx/L). Calculating the force F = —dW/dA n , 
we obtain 



F = 
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If the excess charge density is uniform, p(x) =e/L, F = 
for all modes. Incorporation of interactions that are 
screened at distances longer than the distance to the gate, 
yields a force that scales as L~ 2 . In this case, the coupling 
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parameter g also scales as L" 2 and typical values are in 
the order of 10 -3 , in apparent contradiction with the 
experimental data. 

Assuming that the charge is localized in the center 
of the tube, p(x) = eS(x — L/2), the force is zero for 
odd harmonics, but for even harmonics, n = 2Z, it reads 
Fi = (—l) l e 2 po27rl/L. This results in a coupling param- 
eter that is length independent and scales as higher 
modes are coupled weaker to electrons. Numerical esti- 
mates show that g ~ 1. Localization of an electron in a 
point away from the middle produces coupling to both 
odd and even modes. Note that the electron does not 
have to be strongly localized to produce a g ~ 1. Such 
a non-uniform density can be created by impurities lo- 
cated in the substrate, or induced by a redistribution of 
electrons in a suspended tube bent by a underlying gate 
electrode 30]. 

An interesting feature of the data is the appearance of 
negative differential conductance (NDC) in the current- 
voltage characteristics. NDC is very pronounced in 
Fig. 3a, but also present in Fig. 3b, c. Although several 
explanations for NDC have been put forward, its origin 
remains unclear. Koch and von Oppen ^} showed that 
for low relaxation and strong e-ph coupling, NDC fea- 
tures appear, although they do not follow regions with 
strong positive differential conductance (PDC) as in our 
data. McCarthy et al. p have shown that NDC features 
can be due to an e-ph coupling that is voltage depen- 
dent. Their calculations also reproduce the catastrophic 
current decrease of Fig. 3a for bias voltages higher than 
3 mV. However, at the moment we do not know how such 
a voltage dependence would arise in suspended tubes. 
Finally, Nowack and Wegewijs [2^ have considered a 
Franck- Condon model with a coupling to an electronic 
and its excited state. They show that the competition be- 
tween the two states generates strong NDC effects. NDC 
and PDC lines may have the same gate voltage depen- 
dence preceded by a region of suppressed current. This 
scenario may especially be relevant for the data in Fig. 
3a. 

In summary, transport measurements on suspended 
SWNTs show signatures of phonon-assisted tunneling, 
mediated by longitudinal vibrational (stretching) modes. 
The current-voltage characteristics show multiple steps 
whose heights are in reasonable agreement with the 
Franck- Condon predictions if the e-ph coupling constant 
is of order unity. Suspended nanotube quantum dots 
form an interesting model system for future studies on 
the interaction between single electrons and quantized 
phonons in the intermediate to strong electron-phonon 
coupling limit. 
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